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Abstract. We develop an algorithm to compute Fourier expansions of vector val- 
ued modular for Weil representations. As an application, we compute explicit linear 
equivalences of special divisors on modular varieties of orthogonal type. We define 
three families of Hecke operators for Jacobi forms, and analyze the induced action 
on vector valued modular forms. The newspaces attached to one of these families 
are used to give a more memory efficient version of our algorithm. 
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1. Introduction 

Let £ be a lattice of signature (2, n). One can define a modular variety Xc attached 
to C (see Section [9]). The Picard group Pic(V£) of Xc contains a subgroup PiCsp(V£) 
that is spanned by so-called special divisors Z{m,^). In the past decade, we learned 
how to compute the rank of PiCsp(V£) under certain assumptions on C Most recently, 
this problem was resolved by Bruinier |Brul2| (see |Bru02| for earlier results) in the 
case n > 2 assuming that C splits a hyperbolic plane U = ( ? J ) and a rescaled copy 
U{N) of U. In this setting, the rank of PiCsp(V£) essentially equals the dimension of 
a certain space of vector valued modular forms. In this paper, we develop a method 
to compute relations of special divisors in PiCsp(V£). 
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Theorem 1.1. Let C = U ® U{N) © £'(—1), where CJ is positive definite. Then 
there exists an algorithm which computes linear equivalences between any finite set of 
special divisors Z{m,fi). 

For example, for C = (8) a lattice of rank 1, we have 

Z(i|,3) = 2Z(^,l)-2Z(^,3). 
More such relations will be given in Section O 

1.1. Special Picard groups and vector valued modular forms. The computa- 
tion of Picard groups is intrigued, and explicit determination of linear equivalences of 
divisors is even more difficult. In the case of orthogonal varieties Borcherds provides 
a construction of meromorphic functions, called Borcherds products, whose divisors 
are explicit sums of special divisors jBor98] . If n > 2, there is hope that the converse 
of Borcherds result holds. Every meromorphic function whose divisor is a linear com- 
bination of special divisors should be a Borcherds product. As we have mentioned, 
in some cases this converse theorem is proved. 

Borcherds construction works as follows. Given a weakly holomorphic vector valued 
modular form of type pc (the Weil representation of the discriminant form associated 
to C) consider its Fourier expansion. 

a(m, /i) exp(27ri mT)e^. 

At mGQ) 

Here the components are index by /x and the are certain basis elements for the 
representation space of pc- If a{m, /x) G Z for m < and a(0, 0) = 0, then there is a 
meromorphic function on Xc whose divisor equals 

^ ^ a{m,p)Z{m,p). 

At 0>mGQ 

In the light of the converse theorem, we are thus reduced to considering principle 
parts of weakly holomorphic modular forms, when studying PiCsp(X£). 

Serre duality can be used to determine which principle parts occur |Fis87t IBorOO] . 
This brings into focus holomorphic vector valued modular forms of weight k = and 
type Pc = Pc{-i) (the dual of pc), the space of which we will denote by Mk{pc). Since 
£ is indefinite these spaces are difficult to compute. In the literature, one can find two 
methods which allow to determine Fourier expansions of elements of Mk{pc)- First of 
all, there are Eisenstein series, whose coefficients can be calculated efficiently |BK01| . 
If C splits two hyperbolic planes, similarly efficient calculations can be done for theta 
series |Sie51] ( fBor98] gives explanations in terms of vector valued theta series). The 
drawback of relying on these constructions is that Mk{pc) is often not spanned by 
Eisenstein series and theta series. A second method, which was brought to attention 
by Scheithauer [Sch04[ [Sch08] . is to average usual modular forms. One can show that 
every vector valued modular form arises this way. The drawback in this case is that 
one has to explicitly know Fourier expansions at all cusps of modular forms for groups 
of non-squarefree level. This is an issue which is not easy to settle, and which, for 
example, in [ Raul2b] was dealt with by employing eta products. In general, averaging 
does not seem well suited for large scale computations. As a side node, we mention 
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that there is a third method, that apparently is not treated in the hterature. One can 
use Eisenstein series and theta series for harmonic polynomials in order to construct 
elements of Mk^ii2{pc)- Suitable linear combinations will have vanishing order greater 
than or equal to /, and this vanishing order is sufficient to obtain elements of Mk{pc) 
after dividing by A', where A is the discriminant modular form. 

Summarizing, we can say that averaging as described by Scheithauer is the only 
systematic way to span Mk{pc)- This method, however, relies on another difficult 
computational problem. We will therefore use Jacobi forms to compute Mk{pc)- 

1.2. Vector valued modular forms and Jacobi forms. If C was positive definite, 
for every > A; G |Z there would be a well-know link between Mk{pc) and Jfc_,_2+n £, 
the space of Jacobi forms of weig ht k + ^ and index C (see |Zie89] and Section [3] 
for definitions and details). However, C is never positive definite. Section [2] contains 
the proof of the following theorem. 

Theorem 1.2. Let C be an even lattice. Then there is an even positive definite 
lattice CJ such that C®tl\ = C (B U2, where lAi and U2 are unimodular lattices. In 
particular, we have pc = Pc'- 

After applying this theorem we are reduced to the case of vector valued modular 
forms associated to positive definite lattices. For the rest of this section, we let C he 
a positive definite, even lattice of rank N . 

The theta decomposition for Jacobi forms, which is revisited in Section [3l allows us 
to consider Jfc,£, k > ^ instead of vector valued modular forms. Note that Jk,c is zero 
except if k is integral. It is useful to distinguish two types of Jacobi forms. Jacobi 
forms of scalar index, that is those whose Jacobi index is a lattices of rank 1, are by 
now classical. They were studies systematically in |EZ85| . but first instances of such 
forms already appear in Jacobi's work. In Skoruppa's thesis |Sko84] . an algorithm 
to compute such Jacobi forms was given. Together with the newform theory, hinted 
at in |EZ85] and proved in [SZ88, .SZ89J , this provides us with an efficient means to 
calculate Fourier expansions of Jacobi forms, even if their weight and index are large. 

The case of lattices with > 1 was, so far, not considered systematically. For 
applications that appeared in literature, it was sufficient to use special construction, 
which rely on theta series - see for example [CGO8]. In Section HJ we discuss an algo- 
rithm which gives us Fourier expansions of a basis of J^ £ up to arbitrary precisions. 
The idea behind this algorithm is as follows. Restricting a Jacobi form 0(r, z) to 
z = s^^z', 2;' G C, s G we obtain a Jacobi form of the same weight and scalar in- 
dex. We first prove that sufficiently many of these restrictions are enough to uniquely 
determine (p. In a very special case and for some fixed weights and indices, a similar 
result was used by Das |DaslO| to analyze the structure of so-called hermitian Jacobi 
forms. 

As a second step, we prove that one can construct Jacobi forms this way. Note 
that there is an analog of the restriction z = s^^z' for formal Fourier expansions. 
Suppose that (p is formal Fourier expansion that has all symmetries inherent to Fourier 
expansions of Jacobi forms of weight k and index C If sufficiently many restrictions of 
(j) are Fourier expansions of Jacobi forms of correct scalar index, then is the Fourier 
expansion of a Jacobi form. Algorithm 14.31 makes us of this fact to successively 
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compute restrictions, until the only formal Fourier expansions that restrict to Jacobi 
forms are those which come from Jfc_£. 

We have to remark that our Algorithm depends on enumerating short vectors 
in C In order to use our algorithm we only need to enumerate the lattice once 
for all weights, as long as certain conditions are met. Actually, L has rank 6 or 7 
for all computations that we performed. These ranks are feasible even with non- 
parallel implementations of lattice enumeration that are not optimal. In this sense, 
our algorithm is good enough to yield new results even without consideration of 
technical details of the implementation. In more complicated cases, £ can be too large 
to enumerate sufficiently many vectors so that Algorithm l4.3l finds enough restrictions. 
At the end of Section [7] we comment on further developments that, we think, can 
mitigate the impact of large lattice ranks on performance. 

An idea which similar to the one behind our algorithm was employed by Poor and 
Yuen in order to compute Fourier expansions of Siegel modular forms of genus 4 |PY07] . 
Restricting Siegel modular forms to many different elliptic modular curves, embedded 
into the Siegel upper half space, they obtained relations between Fourier coefficients. 
When increasing the number of such restrictions these relations become sufficient to 
determine the space of Fourier expansions that come from Siegel modular forms inside 
the space of all possible Fourier coefficients. Poor and Yuen could not prove that this 
is always the case, and such a proof is not available until today. Their algorithm, 
however, works perfectly fine in practice. 

There is one fundamental difference between our algorithm and the one in |PY07] . 
While Poor and Yuen would have to increase the number of restrictions if they wanted 
to compute Fourier expansions up to higher precisions, this is not the case for Jacobi 
forms. The reason for this is that (13. 7p and (13. 8p hold, while there is no correspond- 
ingly powerful tool in the case of Siegel modular forms. 



1.3. Hecke operators and newform theories. Hecke theory and the resulting 
theory of newspace can be used to speed up computations of classical elliptic modular 
forms. Section [5] contains investigations related to this problem. In [SZ89] . it was 
proved that we have a direct sum decomposition 



of the space of Jacobi forms of index £ = (2m). This decomposition was already used 
in |Sko84| to provide an optimization for the computation of Jacobi forms of scalar 



index. We generalize the Hecke operators Vi and Vi to operators Vs and Vi on Jacobi 
forms of arbitrary indices. The definition is straight forward, but seemingly it is not 
given in the literature. One might hope to find an decomposition analog to the above 
for all Jacobi indices £, but this is not the case. The lattice ( I § ) provides us with a 
counter example in the following way. Assuming that a direct sum decomposition into 
newspaces exists, it is not hard to write down a recursive formula for the dimensions 
of newspaces. Computing these hypothetical dimensions of newspaces, one expects 
to obtain non- negative values, but if A; = 4 one obtains —1. We will consider the 
question which decomposition yields a direct sum in a sequel. 
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We investigate the induced action of both generahzations Us and Vi on vector valued 
modular forms. The action coming from Ug coincides with the action of an operator 
that was already studied by Scheithauer and Bruinier [Schllt [Brul2| . However, the 
induced action of Vi is new. It corresponds to a curious extension of the discriminant 
group, and deserves further investigation. As it is only loosely related to the theme 
of this paper, we defer this study to a future publication. 

The main result of Section E] says that newforms with respect to the Hecke operators 
Us can be characterized by restrictions to Jacobi forms of scalar index. As a result, 
we are able to describe an Algorithm which makes use of the action of Us to reduce 
memory consumption of computations of Jacobi forms. Since at this stage we are 
not able to prove a newspace decomposition for Jacobi forms which is a direct sum, 
this does not result in a reduction of run time. It does, however, enables us to split 
up computations. This is an useful feature when computing large spaces of Jacobi 
forms. 

1.4. A reference implementation and two families of examples. We provide 
an implementation of Algorithm 14.31 in Purple Sage |S"'"lla] . Section [7] contains in- 
structions on how to install and use it. In Section [71 we also discuss shortcomings of 
our implementation and possible improvements. Regardless of these imperfections, 
our implementation can be used for basic computations. It is also meant to be a 
reference implementation, which clarifies Section |H 

The reference implementation allows to compute Fourier expansions of Jacobi forms 
of arbitrary weights > 2 + y. Fourier expansions of any other weight and of weakly 
holomorphic Jacobi forms can be obtained by division by powers of the modular 
discriminant. The Fourier expansion framework which was introduced in [Raull] 
provides facilities to manipulate and perform linear algebra on Fourier expansions. 
We illustrate how to invoke the corresponding commands and give examples. 

In Section [HI we provide Fourier expansions of vector valued elliptic modular forms 
for cyclic discriminant forms of even order. If these discriminant forms come from 
positive definite rank 1 lattices, the attached vector valued modular forms can be 
easily computed using Jacobi forms of scalar index. If the associated lattice of rank 1 
is negative definite, then they correspond to skew-holomorphic Jacobi forms, whose 
Fourier expansions are not easily accessible. In this case, only our new algorithm 
allows to perform explicit computations. Section [8] contains some printed tables of 
Fourier expansions. Fourier expansions for weights between | and ^ can be found 
at the author's homepage |Raul2a] . 

1.5. Further Applications. The application to special Picard groups that was ex- 
plained at the beginning of the introduction can be found in Section [9l At this 
place, we briefly discuss two further applications, that we will not treat in detail. 
The first concerns explicit computations of Borcherds products, that was treated 
in |GKR11] . We have already explained Borcherds's construction in Section 11.11 
This construction does not only work for meromorphic forms, but in general a(0, 0) 
may be an arbitrary integral number. The resulting Borcherds product will have 
weight a(0,0) / 2, and it will be holomorphic if a{n, fi) > for all n < 0. Such 
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Borcherds products are important in several areas, including string theory |JS05| and 
M24-moonshine |EUTlll [UDT2] . 

Recall that we assume that L has signature (2, n), n > 2. Borcherds products arise 
from weakly holomorphic vector valued modular forms of weight and type pc- 
If, for example, £ is a lattice of the form U ®U ® £'(—1), then pc = pc- Since 
C is necessarily positive definite, we can apply the theta decomposition without 
problems. In order to obtain weakly holomorphic Jacobi forms, which serve as input 
data to Borcherds's multiplicative lift, one can compute vector valued modular forms 
of weight + 121, for some < / G Z. Dividing the resulting Fourier expansion by 
A', where A is the discriminant modular form, gives us a basis of weakly holomorphic 
vector valued modular forms with pols at infinity of order less than or equal to I. The 
needed computations of Fourier expansion of Jacobi forms of index C are rather easy 
for, say, n < A. 

The second application of vector valued modular forms which we want to discuss 
is based on reversing the ideas in |Sch08] . Scheithauer pursued a detailed study of 
averages of scalar valued elliptic modular forms. The Fourier expansions of resulting 
vector valued modular forms depend on the Fourier expansions of the input data 
at all cusps. Modular forms that one starts with typically do not have squarefree 
level. As a consequence, it is quite demanding to obtain expansions at cusps different 
from oo. On the other hand, SL2(Z) acts transitively on the set of cusps C of a 
congruence subgroup F C SL2(Z). Denote the associated permutation representation 
by p. A modular form / of weight k for F gives rise to a vector valued modular form 
of weight k and type p. Its components with respect to the obvious basis of C[C] are 
simply given by the expansions of / at the various elements c of C. 

Since p is unitary, it can be decomposed into a direct sum of Weyl representations. 
Our work thus allows us to compute the Fourier expansion of the vector valued 
modular form with prescribed /oo- Thus we obtain all f^. 

Acknowledgement: The author thanks Max Plank Institute for Mathematics, Bonn, 
Germany for granting access to its computation server. 

2. Discriminant forms and positive definite lattices 

In this section, we relate discriminant forms and positive definite lattices. In par- 
ticular, we explain Algorithm 12. 3[ which allows us to find a positive definite lattice 
that represents a given discriminant form. A useful overview over discriminant forms 
can be found in |Nik79| . 

Definition 2.1. Let D be a finite abelian group that is equipped with a quadratic form 
£)^Q/Z, 7h-)'7^/2. We call D a discriminant form if the induced bilinear form 
(7, 6) := (7 + 5)^ — 7^ — 5^ (mod 1) is non- degenerate. 

In |Wal63] . is was shown that all discriminant forms arise from lattices. Even more 
was shown: Two even, non-degenerate lattices are stably isomorphic if and only if 
the corresponding discriminant forms are isomorphic. To illustrate the corresponding 
construction, fix a non-degenerate, even lattice C C and write (• , ■)c for the 
induced even bilinear form. Set qc{,v) = {v,v)c/ 2. The dual lattice consists of 
all vectors w G £ (g) Q such that (f , w)c G Z for all w E C Then the discrimiant 
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form disc L attached to £ is defined as / C equipped with the quadratic form 
qn (mod 1). For many purposes, in particular in the context of Jacobi forms, it is good 
to identify with Z^. Given a Gram matrix L of £, we find that L then corresponds 
to Lll^ . The quadratic form qn on and C is given by f i— )■ -/^^^[f] = 

We start with a basic lemma on representation numbers of the unimodular lat- 
tice £"8. 

Proposition 2.2. Given any positive, even number 2n, there are vectors v,w & Es 
of length 2n and scalar product {v,w) Es = 2n — 1. 

Proof. In order to prove the lemma, it suffices to show that the [2n-i "^^n^)'^^ coef- 
ficient of the theta series attached to Eg, is non-zero. By |Fre91j . this thet a series is a 
Siegel modular form of degree 2 for the full modular group. By results of Igusa |Igu62| , 
it is a (non-zero) multiple of the Siegel Eisenstein series E^ . The Fourier coefficients 
of El are invariant under unimodular transformations, so that it suffices to analyse 
the Fourier coefficients of 

\-l l) \2n-l 2n J \ ly^l 2) ' 

(2) 

That is, we have to analyze the first Fourier Jacobi coefficient of ii^4 . It is a weight 4, 
index 1 Jacobi form in the sense of |EZ85) . Since E^^ is a Maafi lift, this Fourier 
Jacobi coefficient is non-zero. And since the space of weight 4, index 1 Jacobi forms is 
generated by the Jacobi Eisenstein series E^^i, we are reduced to analyzing the Fourier 
coefficients of this Jacobi form. Its coefficients were studied in |EZ85] . Theorem 2.1. In 
particular, there is a remark preceding the theorem saying that the Fourier coefficients 
of ^4 

1 are positive. This proves the proposition. □ 

The main result of this section is the following algorithm and Theorem 11.21 

Algorithm 2.3. Given a lattice C the following algorithm computes a positive definite 
lattice C that is stabily equivalent to C 

(1) Set C ^ C As long as C is not positve definite iterate the following steps. 

(2) Ghoose a vector v E C of negative length which is as short as possible. Set 
n i qc'iv). 

(3) Ghoose a pair of vectors Wi,W2 G Es satisfying qEg{wi) = qEs{w2) = n and 
{wi,W2)es = 2n-l. 

(4) Set C ^ span(t' + wi, v + ^2)^ C C ® E^. 

Proof of Theorem \1.2[ By what we have observed above, we have D = disc£' for 
some even, non-degenerate lattice C. Applying Algorithm 12.31 we get the desired 
lattice C □ 

Corollary 2.4. Let C be a discrimiant form. Then there is an even non-degenerate 
positive definite lattice C such that D = disc£. 

Proof of correctness and termination of Algorithm \2. ^ Whenever Step [2] is done, £' 
is not positive definite. In particular, there is some vector of negative length in C . 
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By Proposition l2.2[ vectors Wi and W2 as in Step[3]always exist. The lattice span(f + 
^1,^ + ^2) is a hyperbolic plane. Indeed, we have qc's>Es{v + Wi) = qc'{'v) + qEs{wi) = 
1, and qc'(BEs{v + W2) = 0. Further, {v + Wi,v + W2)c'(bEs = 2g£'(w) + (wi, ^2)^8 = -1- 

From this we see that the number of negative eigenvalues of span(f + f + ^2)^ C 
C © Eg is one less than the number of negative eigenvalues of C. This proves that 
Algorithm 12.31 terminates and that the output C is positive definite. □ 

One should remark that the algorithm is not optimal in the sense that C possibly 
contains a unimodular lattice. 



Example 2.5. The lattice C 



/2 
1 

1 
1 

-1 

V-1 



2 1 ^ 
1 2 . 

1 

2 
1 
1 



and 

1 
1 
2 


-1 
-1 



1 
1 

2 




-1 




2 
1 



-1\ 

-1 

-1 
-1 
1 

2/ 



are stabily isomorphic. Since the discriminant group of both lattices is cyclic of prime 
order, we can easily see this by computing the determinant of both matrices. Since 
C has rank 6 it does not split a hyperbolic plane. Algorithm \2. 3[ gives a lattice C" of 
rank 14, which is necessarily of the form C" = C (B Eg. 



3. Vector valued modular forms and Jacobi forms 
Given N e Z^q, let 

Hi,7v := e X 

be the Jacobi upper half plane, where M. = {t = u + iv : f > 0} C C is the Poincare 
upper half plane. It is acted on by the full Jacobi group 

PJ := SL2(Z) X (Z^ ® Z^), 

where the semidirect product is defined via the natural action of SL2(Z) on Z^. We 
write 7'^ = (7, (A, fi)) for a typical element in F'^. Here, A, e are the columns of 
the second component of 7'^, which is viewed as a matrix via the natural identification 
Z2 ® Z^ ^ M7v,2(Z). As before, we write 7 = ^) for a typical element in SL2(Z). 
The action F'^ x Hi tv ^i,n is given by 

'ar + b z + Ar + /i ' 



(3.1) 



^ cr + ' CT + d 

This action extends to a family of actions on function : Hi jv — ^ C. They depend 
on both a weight A; G Z and a (Jacobi) index L G Mi^[Z) which is symmetric and has 
even diagonal entries. In order to define them, denote z^'^Lz by L[z] {z G C^). 

(cr + rf)-^e( -cL[z + \t + ^] /(cr + d) + tL[\] + 2\''Lz) 0(7-^(r, z)) . 
Here and throughout the paper, we use the notation e(x) = exp(27rza;). 
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The following definition, which is folklore, includes a special case of a definition 
made in |Zie89] . where (nonweak) Jacobi forms of higher degree were defined. It also 
generalizes the definition of weak Jacobi forms made in |EZ85] (the index m that was 
used there corresponds to L / 2 in our notation if L has rank 1). 

Definition 3.1. A weakly holomorphic Jacobi form of weight A; G Z and index L is 
a holomorphic function : Hi — ^ C that 

(i) satisfies 0|^^7'^ = for all 'j'^ G F"^, and 

(a) has growth 0(r, ar + (3) = 0(e"("'^^^) for all a,(3 G and some a(a, (3) > 0. 

A weakly holomorphic Jacobi form that satisfies (f){T, ar + /3) = 0(1) for all a, /3 G 
is called a Jacobi form. 

We denote the space of such weakly holomorphic Jacobi forms of weight k and index 
L by j|j ^. The space of Jacobi forms is denoted by Jk,L- 

We now define vector valued elliptic modular forms. The metaplectic cover Mp2(Z) 
of SL2(Zi) is the preimage of SL2(Z) in Mp2(M), the connected double cover of SL2(M). 
Write 7 = ("^) for a typical element of SL2(M). The elements of Mp2(M) can be 
written as (7, r 1— )■ \/cr + d), where the first component is an element of SL2(M) and 
the second is a holomorphic function on H. Since there are two branches of the square 
root, this yields indeed a double cover of SL2(M). 

Given a representation (p, Vp) of Mp2(Zi), /c G |Zi and / : EI — )■ V^, we define 

(/L.„ --M-"p(h.-))-'/(^) 

for all (7,w) G Mp2(M). 

Definition 3.2. Let (p, Vp) be a finite- dimensional representation. A weakly holo- 
morphic vector valued modular form of type p and weight k E ^Ij is a holomorphic 
function / : EI — > such that the following conditions are satisfied: 

(i) For all 7 G Mp2(Z) we have f\k,p1 = f ■ 

(a) We have ||/(t)|| = 0(e"^) for some a > as y ^ 00, where \\ ■ \\ is some 
norm on Vp. 

We will focus on Weil representations attached to discriminant forms D. In this 
case, the second condition translates into a condition on the Fourier expansion. It is 
of the form 

a{n)e{nT), a{n) G Vp. 

-oo<^^e(Q 

In order to define the Weil representation, write 

T:=((il),l) and 5 := ( ( ? ) , v^) 



for the generators of Mp2(Z). The root in the definition of 5* is the principal branch. 
The Weil representation is defined on the group algebra C[D] with canonical basis 
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elements (7 G -D). The images of T and S under are as follows |Sko08] : 

(3.2) p^(r)e^:=e(7')e^, 

(3.3) pz.(5)e,:=— -1^5^e(-(7,5))e,, ^(D) := ' ■ 

We denote the dual Weil representation hj p^,. Given a lattice L, we write p^ and 
PL for pdiscL and pdiscL, respectively. 

Our interest in Jacobi forms stems from the next Theorem. It relates Jacobi forms 
and vector valued elliptic modular forms via theta series. Given / G / LZ^, define 

(3.4) 9LAr,z):= Yl g^'t^'^^V- 

r=l (modLZ^) 

Theorem 3.3 f |Zie89l p. 210]). Let G 4^. Then there is a vector valued weakly 
holomorphic modular form (/i«)igdiscL of type p^ and weight k — such that 

(3.5) 0(r,z)= Y (^L,i{r,z)hi{r). 

igdisc L 

Conversely, given a vector valued weakly holomorphic modular form of type p^ o,nd 
weight k — ^, h3. 5]) defines an element of J'j^^. 

We denote the induced bijection between Jj^^N ^ and yikipb) by 0^, suppressing 
the weight. 

We start revisiting several basic facts about Jacobi forms, which we will need 
to formulate Algorithm 14.31 First, recall the dimension formula for Jacobi forms, 
which follows immediately from |BorOO] , where a dimension formula for vector valued 
modular forms was given. In order to state the next theorem, note that vector valued 
modular forms of weight k — ^ and type pi are supported on the span of + e__^, 
/i G discL, if A; = (mod 2), and on the span of — e_^, if A; = 1 (mod 2). They 
vanish in all other cases, and hence we shall assume that A; G Z. Denote the restriction 
of PL to span{e^ + e_^} or span{e^ - e_^}, respectively, by pk,L- 

Given a unitary matrix M with eigenvalues e(/3j), < /3j < 1, set a{M) = 

Theorem 3.4 ( |BorOO] . Lemma 7.3 and |Fis87] ). //fc > 2 + the dimension of3i:,L 
equals 

Remark 3.5. For reasons of performance, it is not feasible to implement the above 
dimension formula without major modifications using exact arithmetic. The imple- 
mentation of Algorithm \4.3\ discussed in Section Q makes use of the fact that the 
eigenvalues ofpk,L{S) andpk,L{ST) are twelfth roots of unity. Upper bounds for each 
multiplicity of these eigenvalues are computed using interval arithmetic and House- 
holder transformations. Correctness of the computations can then be proved a poste- 
riori by comparing the sum of all multiplicities and the dimension ofpk,L- 

To state the next lemma, which is folklore, set L[s] = s^^Ls for any pair of matrices 
L and s of compatible size. 
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Lemma 3.6. Let cf) he a (weakly holomorphic) Jacobi form of weight k and index L. 
For s G Mjv,iv'(Z) ^A^' G Z^q), the pullback (j)[s]{T,z) = (j){T,sz) (z € C^' ) is a 
(weakly holomorphic) Jacobi form of weight k and index L[s]. 



Proof. The conditions in Definition 13.11 can be verified by a straight forward compu- 
tation. □ 

A weakly holomorphic Jacobi form has a Fourier expansion of the following form: 

(3.6) 0(r,z)= J2 

where q = exp(27rir) and ("^ = exp{2ni YliLi "^i^i)- Notice that the Fourier coefficients 
c(n, r) of a Jacobi form vanish if 2|L|n — L*[r] < 0, where \L\ := det(L) and L* : = 
\L\- is the adjoint of L. 

Because Jacobi forms are invariant under the action of (/2, (A,0)) G F^^, A G Z^, 
I2 the identity matrix of size 2, the Fourier coefficients c(r;,,r) satisfy the following 
relations, which are connected to the statement of Theorem | 

(3.7) c(n, r) = c(n + L[\] / 2 + r''\ r + LA) 

(3.8) c{n,-r) = {-if c{n,r). 
It is convenient to define 

c(A.r):=c(^±^.r). 
By the above relations, c only depends on r (modLZ^). 



4. Computing Jacobi forms 

Because of (13. 7p and ( 13. Sp . the Fourier expansion (13. 6p of a Jacobi form is deter- 
mined by all c(n, r) where r lies in a set TZ of representatives of (Z^ / LZ'^^ /{±1}. 
We can assume that the elements r E TZ have minimal norm among all elements of 
r + LZ^. Set 

neigh(r) = {r' G r + LZ^ : L[r'] = -L[r]} and mult(r) = #neigh(r). 

We call a finite set iS C Z x LZ^ a complete set of restriction vectors (for TZ) if 
the matrix 

F{n,S) := {mult{r)6s^rr=f)^^T^^^,,^^s 

has rank i^TZ. We say that S is minimal if the above matrix is invertible. Obviously, 
every complete set of restriction vectors contains a minimal one. By abuse of notation 
we will write S for the projection of S to its second component. 

Lemma 4.1. Given L and a set TZ C Z^, there is a complete set of restriction vectors 
for IZ. 



12 



MARTIN RAUM 



Proof. It suffices to prove that for each r & TZ there is (f, s) G Z x LZ^ such that 
gtr^ _ ^ g^j^j J, jg ^YiQ unique element in TZ with this property. Choose s' G whose 
components are linearly independent over Q. Set 

e = min {\t^i — V' l)i 
where tr' = s'^^t' , r' G 7^. Choose an approximation s" G Q'^ of s' such that 

<e/(^ max (|r;.|)) 
i<i<-'v 

for all 1 < 2 < A^. Then all s"r' E Q, r' E TZ are distinct, and we can set s = 
s" ■ Icmjdenominator of s^'}i<i<Ar. □ 

For later use, we need strong sets of restriction vectors S^^'^^r). 

Lemma 4.2. Let r E TZ U —TZ. Then there is a finite set iS'^*''(r) that spans Z^ 
and has the following property. Whenever s^^r = s*^r' for some s E S{r) and some 
r' G IJ^g^ neigh(r) U neigh(— r), then r = r'. 

Proof. The proof of lemma I^TTl gives us a vector s which has the second property. Fix 
a basis Vi, . . . ,Vn of Z^. Fix some h > maxi^ |f*'^(r — r')|, where r' runs through 
|J,^g^ neigh(r) Uneigh(— r). Then we can choose 

S^'^'ir) := {s, hs + vi,...,hs + Vn}- 

This set obviously spans Z^, since {vi, . . . , f^v} does. Since 

\{hs + Vif{r' - r)| > h\s''{r' - r)\ - (r' - r)| > 0, 

its elements satisfy the second property given in the lemma. □ 

We are now in position to formulate an algorithm that computes the Fourier coef- 
ficients c(n, r), n < -B G Z>o associated to Jacobi forms of index L G Mjv(Z), N >2. 
To give a rigorous description, set 

I{L,B) = {(n,r) GZxZ^:0<n<5, 2\L\n - L*[r] > O}, 

and let J^£l,b ■ 0^ h,L ^ C^(^'^) be the map that sends a Jacobi form to its 
truncated Fourier expansion. We denote the subspace of C^^^'^') whose element satisfy 
the relations O and ED by C^iL,By,r^ ,k ^ Given v G C^(^-^), let v[s] G C^(^H-S) be 
the vector, whose (r;,,f)-th component is given by Ylrez'^ s^^r=f'^{n,r)- This formal 
pullback map is made so that = J^£Lt^s\,B{'P[s\) for any cj) G Jfc,L; A; G Z. 

Algorithm 4.3. Fix k E Z, L as above, and < i? G Z. The following algorithm 

computes FSL,B{^k,L) ■ 

1: Compute a set 7Z of representatives r of (Z^ / LZ^) / {±1} with minimal norm L[r] 
in r + UL^ . 

2: Compute a minimal complete set S of restriction vectors for TZ 
3: Set B ^ max {B, max^{dim(Mfc+L[s], dimMk-i+L[s]}} ■ 
4: For each s eS, compute Jk,L[s\,B ■= J^SL[s],B{h,Ll4) ^ Q^'^'^^''''^^- 
5: Compute Jk,L,B := {{v[s]),es : v G Q^iL,By,r\k-^ ^ 0,g5Q^(^W'^^ 
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6: If the dimension of Jk,L,B H ^^^g Jk,L[s\,B does not equal div[iik,L, enlarge S by 
the shortest possible vector with respect to qc ■ Continue with Step O 

7; Compute the inverse image of Jk,L,B^@ s&s •^k,L[s\,B under v ^ {y[s\)si^s- It equals 
TEb ( J^.l) ■ Truncate it to the required precision. 

Remarks 4.4. 1) An implementation of this algorithm in Purple Sage |S"'"lla 



zs 

available at |Raul2c] . For reasons of speed, this implementation represents Jk,L,B^ 
®seS'^k,L[s\.B o-s the intersection of three spaces. The reader is referred to the 
implementation for details. 
2) A precise set S which is necessary for Step\^to pass, can be deduced from the proof 
of termination. We do not make this explicit, since in practice any sensible choice 
of S is sufficient. 

Proof of Theorem \l.l\ As explained in the introduction, by the resuhs in [Bor98] 
and |Brul2| . it is sufficient to compute Fourier expansions of vector valued modular 
forms of weight and type pc- Since only finitely many divisors Z{m,fi) are 
considered, these Fourier expansion need to be computed only up to a finite precision. 
Such expansions correspond to Fourier expansions of Jacobi form via Corollary 12.41 
and Theorem 13.31 Now Algorithm 14.31 can be used to compute the desired linear 
equivalences. □ 

Proof of correctness of Algorithm \4.3[ We first prove correctness of Algorithm 14.31 
The Fourier expansion of Jacobi forms in Jk,L is determined by c(n, r), where r 
Since S is contains a complete set of restriction vectors for TZ, the restriction map 

is injective. The condition B > msiXs{dimMk-^L[s]y dimMk-i-^-L[s]} guarantees that 
the J^Sl[s]^b, s G iS are injective. 

We exploit the commutative diagram 

Jk,L ^ Jk,L,B 

®(-)H ®{-)H 

Jfc,L[s] ^ ©se<S J'k,L[s] 

By what we have show the composition of the left arrow with the bottom arrow is 
injective. Hence so is the top arrow. By the dimension check in Step [6], it is an 
isomorphism. This proves correctness. □ 



The proof that Algorithm 14.31 terminates will occupy the rest of this section. We 
will denote the components of z hj Zi, . . . , zn- We say that a function — )■ C is 
elliptic with respect to L and r G EI if is invariant under I (A, /i) for all A, /i G Z 



N 



Here, 

0|^^^(A,/x) := ((r',^)^0(z))|, J/2,(A,/x)) 

for some k G Z — the slash action on the right hand side is independent of weight k. 
We will henceforth use the language of lattices, that was introduced in Section [2l 
Then L will be a fixed Gram matrix of a lattice C. 
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Lemma 4.5. Let (j) he a holomorphic function that is elliptic with respect to a lattice C 
and r G H. Then there is a positive integer b that only depends on L such that if 

(4.1) C -> C"-(ei,^_i), z, ^ 0(r, (^1, . . . , z^)) 

has more than b zeros counted with multiplicity, then = 0. 

Proof. Choose a vector / = {li, 0, . . . , 0)*'' G L^^Z^ fl Z^. Since is elhptic with re- 
spect to C and r, it is, in particular, invariant under the transformations |^ ^(/2, (0, ei)) 
and |^^(/2, (^,0)), where Ci is the first canonical basis vector of Z^. In other words, 
the function fl4.ll) is elliptic in the usual sense. Then the result follows by applying 
the Residue Theorem in the same way as in |EZ85) . Theorem 1.2. □ 

Lemma 4.6. Let cj) he a holomorphic function that is elliptic with respect to a lattice 
C and r G H. Suppose that N > 2. Then there is a set S of matrices s G Mn^n^i^Z) 
that only depends on C such that if 4>[s] = for all s E S, then vanishes identically. 

The following corollary can be deduced from Lemma 14.61 using induction. 

Corollary 4.7. Let he an elliptic function for a finite index suhlattice L' of L, and 
suppose that N > 2. There is a set S of vectors in that only depends on L' such 

that «/0[s] = for all s E S, then = 0. 

Proof of Lemma \4.6\ We will show that, given Zi E Q + iQ, there is a set S'{zi) of 
matrices s G M7v,Ar-i(Z) such that if 0[s] = for all s G S'{zi), then 

0^1 {Z2, ...,Zn) := 0(^1, Z2,..., Zn) = 0. 



Suppose that we have already proved that S'^Zi) exits for all Zi. Apply Lemma [4.51 
to 0, C and r to obtain b. Then we can set S := IJ'^'(^i) where the union runs over 
b distinct ZiS. 

We are hence reduced to finding S'{'zi). In the case of i^i = 0, we can choose 



(0(^-1))^ 



where 0*^^"^^ is the vector of length — 1 whose entries are all 0. To prove the other 
cases, note that 0^-^ is elliptic for some finite index sublattice of C We use induction 
on > 2. We consider the case N = 2 separately. We employ Lemma 14.51 to get 
that can be applied to 0^^. We can set S'(zi) = {s'i} ^^-^i^^- C Q, where the 
are pairwise not equivalent modulo 1. To see that that this choice of S'{zi) works, 
write s'i = Pi / qi with coprime integers pi and gj. We have 

<P[{(l^,PiY']{zi I Qi) = (t){zi,zis[) = (j)z,{zis[). 

Suppose that the left hand side vanishes for all i. Then so does the right hand side. 
By Lemma 14. 5[ we then find that 0j^ = 0. This completes the proof in the case 
of A^ = 2. 

In the case of A^ > 2 we can assume that the statement of the lemma holds for all 
A^' < A^. Given ii, choose S for 0^^ as in the lemma and call it iS(ii). A straight 
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forward computation shows that we can set 

1 (o^^-^))*"^ 



0(^-1) s' J JseSil^) 

□ 



Proposition 4.8. Let (j){T, z) he a function that is 1-periodic m r G EI and an elliptic 
function in z with respect to C and r for all r. Fix a weight A; G Z, and suppose that 
N > 2. Then there is a set S C Z^^^ that only depends on C such that if (f)[s] is 
a Jacobi form of weight k and index L[s] for all s ^ S, then (f) is a Jacobi form of 
weight k and index L. 

Proof. Set S = (I^q)- Using the group structure of F'^, we see that (0 — (p^^^S) 
is elhptic in z with respect to C and r for all r. Choose S as in Corollary 14.71 and 
assume that 0[s] — |^ = for all s G 5. For z' G C^^"*^, we have 

(0 - <Pl^^S) [s]{t, z') = 0[.](r, z') - r-'<p{—, -) = {<p[s] - <l>[s]l^^S) (r, z') = 0. 
Hence, by Corollary 14. 7^ we have — 0| 5 = 0. Since F'^ is generated by S, T, and 



\k,L 

1? ® C F'', this proves that is a Jacobi form of weight k and index L. We have 
thus proved that the proposition holds with S as in Corollary 14.71 □ 

Proof of termination of Algorithm \4.'J\ Given B' < B, let Ub\b = span (e^^^ • B' < 
n < 5) C Q^(^'^), where en,r denotes the canonical basis vector attached to {n,r) G 
I{L,B). Set 



Jfc,i:=lim n j7fe,L,B'+WiJ',B C0^(^'°°) 

0<B I ' 



0<B'<B 



Since dim J'k^i.B' < oo for all B and dim J'k^L,B' < dim J'k.L.B for B' < B sufficiently 
large, it suffices to show that 



J^^oo '■ Jk,L ^ 



k,L 



is an isomorphism. 

The Fourier coefficients of Jk.L[s\ have at most polynomial growth in n, and so do 
the coefficients of any element of Jk,L- Consequently, any such corresponds to an 
absolutely convergent Fourier expansion on Hi tv- By construction these functions are 
elliptic in z and 1-periodic on r. Lemma [4.81 allows us to find a set S^^L) so that it is 
sufficient to check that 0[s] is modular. Since in Step[6]iS is enlarged by short vectors, 
it will eventually contain Sk,L- This proves that Algorithm 14.31 terminates. □ 



5. HECKE OPERATORS AND NEWFORMS 

We define three types of Hecke operators that generalize those studied in |EZ85| 
in the case = 1. The aim of this section is to relate two of them to an action on 
vector valued modular forms. We give the third one only for completeness. The first 
operator is defined for all s G MAr(Z) fl GLAf(Q). The second and third are defined 
for all < / G Z that are coprime to \L\. 
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Expressions for them can be obtained from the parabohc Hecke algebra. In order 
to define it, consider the group 



GSp 



N+l 



M^'JM = iy{g)J, v{g) E Z}, 



{(?GM2(^+i)(Z) : 

^q(N+1,N+1) 

where O^^^'^^) ^ M,v..w.(Z) is the zero matrix. The Jacobi group T'j^ can be embedded 



J 



q{N+1,N+1) 



into Spjv+i(Z) 



{g G GSpjv+i(Z) : iy{g) 



( 



■ a b' 
. c d 1 



(A,/i) 



a 
X 

c 



= 1}: 

0(1.^) 
In 



d 



q(N,N) q(JV,1) 



q{N,N) 

In 



C/i 



It sits inside the monoid GT'j^ of all matrices 



0(1.^) 



where a',d' e M7v(Z)nGL 



a 

X a' fi 

c 0(1-^) d 

Ny'Q), A G Z^, and 



q{N,N) 

d' 



G GSp^+i(Z), 



a b' 

■ c d I 



G M2(Z)nGL2(Q). The action 



of such a matrix, call it g for the time being, on Jacobi forms is 



-(N+\)k 



(cr + d)-^e(^- (cr + d)-^L[z + Ar + /i] + 2X''Lz + L[A]r)0(^^^^, 

The Jacobi Hecke algebra "H^ is the algebra attached to the Hecke pair {T\j, ^^l?} 
- see [KriQOj for details. Note that it is a noncommutative algebra. The operator U'^ 
is induced by the double coset that is generated by the block diagonal matrix 

diag( det(s), det(s)s, det(s), s*) . 

It maps Jacobi forms of index L to Jacobi forms of index L[s\. Given / as above, the 
operator V/ is induced by the set of double cosets which come from the matrices 

/ a 0(1'^) h 0(1'^) \ 

c 0(1'^) d 0(1'^) 
yoWD o(^'^) O(^'i) In ) 

where ( " ^) runs through all matrices with determinant /. This operator maps Jacobi 
forms of index L to Jacobi forms of index IL. The operator T[ preserves the index. 
It is induced by the set of double cosets which come from matrices 

/ a 0(1'^) h 0(1'^) \ 
A Un O(^'i) 0(^'^) 
c 0(^'^) d 0(^'^) 



'AT 



where f ^ d ) I'uns through all matrices with determinant such that gcd 7 is a square. 
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We will use rescaled versions of Ug, V/, and T/. Their actions are defined as follows. 

(5.1) <j)l^^Us = <j){T,sz) 

(5.2) 0L,,V^ = /^ E + 



Iz 



7er\M2(z) 

det ■y=l 



CT + d 



(5^3) ^l^^^r, J: 

det j=fi 

gcd 7=n 

Hecke operators for classical elliptic modular forms will occur later. By abuse of 
notation, we denote them also by T;. 

If is a Jacobi form of index L, then (j)\k,L Us has index L' := L[s]. In the language 
of lattices, the relation of indices corresponds to the inclusion C := sC ^ C On 
the other hand, inclusions of lattices manifest themselves as isotropic subgroups H 
of D' = disc C. We now show that the action of Us on the components of the theta 
decomposition is the same as the action of defined in |Schll] . Recall the definition 
of this operator. Let H G D' he an isotropic subgroup of a discriminant form D'. 
Then, by definition, we have 

for every vector valued modular form / of type p^, with D = H-^ / H. In the present 
situation, we can choose the isotropic subgroup H := C / C G disc£'. We then have 
H^/H = disc£ C disc£'. 

Proposition 5.1. The following diagram commutes: 

T , T 



e 



L 



e, 



Proof. From (15. ip and the theta decomposition of G Jk,L 



we see that it suffices to prove the equality 



(5.4) 



9l,^{t,Sz)= E dL>,s^r^+s^rL^,{T,z). 
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By passing to the Fourier expansion of Ol,^-, this becomes an easy calculation. 



Es*"'At+s"^L/x' (modL[s]Z^) 



□ 



We next treat the operator Vi. In the case = 1, the action of Vi on the Taylor 
expansion of (f) is given by usual Hecke operators. This result can be extended to 
arbitrary N . In order to make a precise statement, we recall the operator given 
in [EZ85] for < G Z. On Jacobi forms cj) of even weight k and index m it equals 



^ \- (-l)A' (2Z/ + 1)! (fc + 21/ - /i - 1)! , ,2.+l-2M, dr 



If k is odd, set 

fi! (2;/ + 1 - 2f,)! (*: + I. - 1)! ^2«' '2m 

In all cases that are not covered by the above definition we have 'D,y{(f)) := 0. 
Proposition 5.2. For all MAr,Ar/(Z), 1 < A^' < A^, u;e have 

In particular, the following diagram commutes for all s G and all < u E 1^: 



Mfc+, Mfc+, 



Proof. The first claim follows immediatelly from the definition of VJ. The second 
claim is then a consequence of the corollary on page 45 of |EZ85] . □ 

The action of Vi on vector valued modular forms is involved. Let C be a lattice 
with Gram matrix IL, that is, a rescaled version of C. When identifying both 
and (jC')* with Z^, then the multiplication by / map induces an embedding of the 
module D = disc C into D' = disc The quadratic form is not preserved. Indeed, 
when identifying fi G disc£ with a corresponding vector in Z^, then //i is unique up 
ILZ^ . Hence In corresponds to a well-defined element of disc£'. We have 

There are further multiplication maps from D to D' , which are only well-defined up to 
submodules of Z^ / /Z*^. Let d\l, then d^ gives an element of D' that is well-defined 
up to elements jj,' in the kernel ker( ■ ^ of the multiplication by d map on D' . 
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Define an operator that restricts Fourier expansions on H: 

(5.5) ^a(n)g"| res((i, n') = ^ a(n) g". 

nGQ nedZ+n' 

With this restriction operator at hand, we can define a new Hecke operator for vector 
valued modular forms. Let D be a discriminant form and D = disc C for some positive 
definite lattice C. Then we set 

(5.6) f\sc,{C) = J2C^)''' E {U\Md,T^D'{jif^ + f^'))){lT/d')u^^^,. 

d 1 1 /iGdisc C 

/i'Gker{ ■d\i^i) 

At this stage, sci{C) does not only depend on D but also on a lift of D to a lattice C 
Proposition 5.3. The following diagram computes: 

Jk,L ^ Jfc,L' 



e 



L 



Br 



SC;(£) 



Proof. We choose (q^), O'd = I, b (mode?) as representatives for 7 G r\MAr(Z), 
det(7) = /. Then it suffices to prove for each d \ I that 



E ( E e(^(L-Mr]/2 + A)))5(A,/.)g^(--^M/^+^)d 

r=ix (mod LZ^ ) b (mod d) 
0<Ae(Q 

= E E £(A,/.)g('^)-^M/^^^\^ 

ti'&Z'^ /dZN r=.L{^+Lti') (modiLZ^) 
A6dZ-L-i[^t+LAt'] /2 

The Fourier coefficients c of a Jacobi forms were introduced in Section [3l The neces- 
sary calculations are analogous to those in the proof of Proposition 15.11 □ 

We end this section with a decomposition into newspaces which is finer than the 
one employed in the next section. Set 

Jfcl°''^= E E Jfc.^'L,L'/i^^^'' 

seMjv(Z) l\L' 
L=L'[s] 1^1 or det(s)5>^l 

Jfc,L - (J/fc,L ) ^ Jfc,L5 

where l\L' means that V / 1 & M7v(Q) has integral entries and even diagonal en- 
tries, and s = In, if det(s) = 1. The decomposition which generalized the newform 
decomposition in |EZ85| in the most straight forward way is 

(5-7) Jk,L = E E 

SGMiv(Z) l\L' 
L=L'\s] 
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where again s = In, if det(s) = 1. We show that this decomposition is not a direct 
sum by giving a counter example. Suppose it was. Then one can recursively compute 
the dimension of i"'^™- Thus recursively define 

sgMjv(Z) l\L' 
L=L'[s] l+\ or det{s)7^1 

where once more s = In, if det(S) = 1. 

Proposition 5.4. Let L = (Ij) and k = 4. Then we have d^J""""^ = -1 < 0. In 
particular, for these values of L and k, /i5. ?)) is not a direct sum. 

Proof. The corresponding code can be found at |Raul2aj . □ 

6. Modification of Algorithm 14.31 

In this section, we give an improved version of Algorithm 14. 3^ which makes use of 
a decomposition of Jk,L- In total analogy to Section 3.1 of |Brul2j, we set 

jfc?= ^>'Mk,L'^s^ 

s6Mjv(Z), dct{s)7^1 
L=L'[s] 

jf/new / jf/oldN^ ^ t 

The orthogonal complement is taken with respect to the Petersson scalar product 

dx dy 



JSL2Z\H ,.^r, 



The central result of this sections says that oldforms can be recognized by their 
restrictions to Jacobi forms of scalar index. 

Theorem 6.1. Let s G M^C^) and assume that L = L'[s] for some V G Miv(Z) with 
even diagonal entries. Then there is an explicit set S C such that the following 
holds. 

For s E S, let Ig be the greatest common divisor of the entries of ss. Suppose that 
7^ G Jfe,L. If(j)[s] G h,L[s]/ii\k,L[s]/i^ s eS, then (p G Jfc.L'l^^^, Us. 



The first assumption of Theorem 16.11 can also be checked on restrictions. 

Proposition 6.2. Let s G M7v(Zi). Then there is an explicit set S C such that 
the following holds. 

For s E S, let Ig be the greatest common divisor of the entries of ss. If 21^ \ L[s] 
for all s E S, then L = L'[s] for some L' G M7v(Z) with even diagonal. 

Proof There is < / G Z such that IZ^ C 5Z^. We show that we can choose 

S = {ir^Ci : 1 < z < iV} U {ir\ei + Ci+i) : l<i< N -1}, 

where Cj denotes the i-th unit vector. We have Ig = P for all s E S, 

Suppose that 21^ \ L[s] for all s E S. We want to show that i^fs""'^] is integral 
with even diagonal entries. Since L[ls~^s'] = PL[s^^][s'] for all s' G Z^, this follows 
immediately. □ 
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Proof of Theorem \6.1[ Recall the definition of TZ and 5'^*''(r), r from Section HI 
We set 

S={j {irh' : s' eS''\r)]. 

Since vectors in S{r) are primitive by definition, we have Ig = I for all / G S. 

Suppose that G •^k,L[s]/i'i\j^ ^gyl■2Ul^ for all s G 5. We have to show that 
4> £ Jfc.L' 1^ ^, f^s- This holds if and only if J-'£((f)) is supported on pairs (n, r) for 
which (s*'')~^r is integral. Recall the Fourier coefficients c and c in Section [31 Because 
L = s^^{L's) and hence c(A,r) = c(A,r + s^^{L'sX)) for A G and also because 
of (13. Sp . it suffices to show this for a system of representatives of (Z^ / LZ^) /{±1}. 
By definition 7^ is such a system of representatives. Henceforth we assume that r eTZ. 

We prove the claim by induction on n. Suppose that all {n',r) with n' < n and 
c{n', r) 7^ satisfy {s~^y^r G Z^. Fix r with c(n, r) 7^ and set f = s^^r. The Fourier 
coefficient of index (n, f) of equals 

^ c(n,r'). 

Using the induction hypothesis, (13.71) . and (13.81) . we can eliminate all contributions 
of c(n',r'), n' < n, r' & TZ. We are then left with the expression 

^ c(n,r). 

r"e-R 

r'gneigh(r") U neigh(— r") 

The sum only runs over r' = r, as can be seen directly from the definition of S^^^{r), 
and hence it is non-zero. By assumption, we have ls\r- We conclude that / 1 s^'^r = 
ls"'{s-^yW. Consequently s"'{s-^y'r G Z for all s' G Since spans Z^, 

we deduce that all entries of {s~^y^r are integral. This proves the theorem. □ 

We formulate an algorithm that makes use of Theorem 16.11 in order to restrict the 
necessary computation to newspaces. 

Algorithm 6.3. Fix k & Z, L as above, and < -B G Z. The following algorithm 
computes J^^L,ij(Jfc,r™) ■ 

1: Compute a set IZ of representatives r of / LZ'^) /{il} with minimal norm L[r] 
in r + UL^ . 

2: Compute a minimal complete setS of restriction vectors forTZ. SetS ^ SVJS'VJS" , 
where S' and S" are the sets given in Proposition ^.^ and Theorem \6.1[ 

3: Set B ^ max {B, maXs{dim(Mfc+2.w, dimMfe_i+L[s]}}- 

4: For each s G S, compute := ^£:lw,s(J^,1h) ^ Q^^^'^''""^ ^Sh,b ■ = 

5: Compute dim J^^^" = dim ^k,L — dim J^^I"^ by recursively applying this algorithm. 
6: Compute J^^^^b := {(^H).e5 : v G Q^(^'^)} C 0,^^ Q^^^H'^). 
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7; Set 

<^k,L'^ ^ "-^^^L ^ ( '^k,L[7]/ri,B I ' 
[ls)s seS' 

where the outer sum runs over tuples {ls)ses of positive integers which satisfy 
gcd{{ls)ses) = 1 and \ L[s] for all s eS. //dim Jj^l'''" ^ dim J^'^*^^, enlarge S 
by the shortest possible vector with respect to qc- Continue with Stepl^ 

8: Compute the inverse image of J'^'^l^^ under v i— )■ It equals J^£B{JkY'") ■ 

Truncate it to the required precision. 

Proof of correctness and termination. Correctness of Algorithm 16 . 3 1 follows from The- 
orem [6?T1 Because we have 

Algorithm 16.31 terminates. □ 

We finish this section with some comments on Algorithm 16.31 First of all, it allows 
for a significant reduction of memory consumption, if the J'k,L[s],B is large, because it 
allows us to split these spaces. In contrast to algorithms for elliptic modular forms 
that make use of the newform decomposition, it does not necessarily reduce the time 
consumption for the following reason. While in the case of elliptic modular forms one 
can show that the corresponding decomposition into newspaces is a direct sum, this 
is not the case for the decomposition of Jk,L arising from the Us- As we have seen, not 
even formula (15. 7p corresponds to such a direct sum decomposition if iV 7^ 1. If an 
analog decomposition could be found, we would be able to compute the corresponding 
dimensions of newspaces recursively. Thus we would not have to compute oldspaces 
in the course of Algorithm 16. 3^ saving, indeed, much time. 

7. Implementing Algorithm 14.31 

This section contains instructions on how to use the implementation of Algo- 
rithm 14.31 that is provided by the author. We will also discuss the algorithm's and 
implementation's bottleneck. 

Install Sage |S^llb] in such a way that you can modify it. Using git |T+12] . clone 
the author's Purple Sage repository at |Raul2c| . 

git clone git : //github . com/martinra/psage . git 
cd psage 

git fetch origin paper_coniput ing_ j acobi_f orms : j acobi 
git checkout jacobi 
sage -sh 
. /build_ext 

Link the downloaded psage library to the site packages folder of Sage's Python. If, 
for example, both Sage and Purple Sage are in the same folder, go there and type: 

PSAGE= ' readlink -f psage/psage' 

In -s $PSAGE sage/local/lib/python/site -packages / 
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At the current stage, it is important that you apply patch 7^13531 at |Rau| to 
Sage. It improves enumeration of short vectors in lattices. Without this patch the 
examples in this section might not even run correctly. You can patch your copy of 
Sage as follows. Go to the directory where you have installed Sage. We assume that 
the patch can be found there, too. 

cd sage/devel/sage/ 

sage -hg patch ../../../ trac_13531 - short_vectors . patch 
sage -b 

Start Sage, either in the terminal or as a notebook, and type: 

from psage . modf orm . j acobif orms import * 
k = 10 

L = QuadraticForm (matrix (2 , [2,1,1,2])) 
B = 3 

jfs = JacobiFormsDl (QQ , JacobiFormDl_Gamma (k , L), B) 

This computes the module of Jacobi forms of weight k = 10 and index L = (12)) 
where the underlying Fourier expansions ^„rC(n, r) are truncated at n < B. 

In order to access the Fourier coefficient at n = 2 and r = (1, 1) of the first basis 
element {Jk,L has dimension 2), type: 

ch = JacobiFormDlWeightCharacter (k , L . matrix (). nrows () ) 

n = 2; r = (1, 1) 

phi = jf s . . f ourier_expansion 

phi [(ch , (n, r) )] 

Currently, multiplication by modular forms is not implemented properly. But by 
using power series instead one can circumvent this short come. We start by extracting 
the Fourier coefficients for all r G 7^, which was described in first paragraph of 
Section HI 

r_reps = ( jfs . f ourier_expansioii_ambieiit (). action () . 

_r_reduced_representatives ) 
R = PowerSeriesRing (QQ , 'q') 
prec = phi . precision 

r_cmps = dict( (r, R ( [phi [ ( ch , (n , r))] if (n, r) in prec else 

for n in range(B)]). 
add_bigoh (B) ) 
for r in r_reps ) 

In order to compute the product of a Jacobi form with a modular form, it suffices to 
take the product of these components with that modular form. 

mf = ModularForms ( 1 , 14) . basis () [0] . qexp (B) 

r_cmps = dict((r, mf * c) for (r, c) in r_cmps . iteritems ( ) ) 

Finally, we have to construct a Fourier expansion from these components. Since we 
want to express the product in terms of standard generators, we can use the attached 
module of Fourier expansions. 
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jfs_phi2 = JacobiFormsDl (QQ , JacobiFormDl_Gamma (k+14 , L) , B) 
phi2_di = {ch:dict( ((n, r) , a) for (r, c) in r_cmps . iteritems () 

for (n, a) in c . diet (). iteritems () )} 
phi2_fe = j f s_phi2 . f our ier _expans ion_ambient ( ) ( phi2_di ) 
phi2 = jfs_phi2( j f s_phi2 . graded_submodule ( None ) ( phi2_f e ) ) 

The result is the following Jacobi form: 

024,l[O] - 24 * 024,l[1] - 255762 * 024,l[2], 

where the 024,l [i] sue the basis elements of J24,l that are computed by our implemen- 
tation. 

Since our major goal is to compute Fourier expansions of vector valued modular 
forms, we illustrate how to convert (p to such an expansion. We start by initializing 
the corresponding discriminant form. 

from psage . modf orm . vector_valued import * 
D = DiscriminantForm (L) 

D . positive_quadratic_f orm == L.matrixO 

The result of this check is True. It indicates that the positive definite quadratic form 
(that is, lattice) corresponding to £ is £ itself. This will always be the case if C is 
positive definite. In all other cases. Algorithm 12.31 should be applied to find such a 
lattice. 

Using the discriminant form, we can extract the Fourier expansion of the vector 
valued modular form which corresponds to 0. 

f = jf s . . theta_decomposition () 

f_simple = dict( ( la . element ary_represent at i on () , c) 

for (la, c) in f . iteritems () ) 

The first command computes the desired Fourier expansion. It will be represented by 
a disctionary whose keys are elements of D. Even though this is good from the point 
of few of abstraction, it is unpractical for storage of such Fourier expansions. It would 
introduce an unnecessary dependence on Purple Sage. Therefor, our implementation 
provides a method that converts elements of discriminant forms to tuples. These 
tuples correspond to vectors in Z^, where D comes form C C Z^. Its entries are 
positive and it is reduced with respect to the row space of L. For example, the 
reduced tuples with respect to {H), whose echelon form is ( J are (0,0), (0, 1), 
and (0,2). 

7.1. Bottlenecks, parallelisation, and extensions. There are two major bottle- 
necks of our implementation that can both be mitigated by parallelisation. When 
computing Jacobi forms for large lattices of small discriminant the most time con- 
suming step is to enumerate short vectors in the lattice. Sage calls Pari |Par03] 
via pexpect, and then sorts the corresponding vectors by length. Obviously, this is 
very slow. It also provides unnecessary data, the vectors' lengths, in some cases. A 
more efficient implementation of Algorithm 14.31 should be based on optimized algo- 
rithms for this task. In particular, it is possible to parallelize enumeration of short 
vectors [DSTO] . 
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When computing Jacobi forms for lattices of large discriminant up to high preci- 
sions, linear algebra can become a bottleneck. Optimization and parallelization of 
this step is of general interest, but it is also possible to split up the original prob- 
lem by using newform theories as indicated in Section [6l We have already made 
clear, however, that currently there no completely satisfactory theory of this kind. In 
particular, it would be important to have explicit dimension formulas for spaces of 
newforms. 

The provided implementation of Algorithm 14.31 does not expose one important 
feature of our algorithm. Since restricting formal Foourier expansions of Jacobi indices 
L to L[s\ for some s G is independent of the weight modulo 2, it suffices to 
calculate the corresponding matrix once. Note that the dependence modulo 2 can 
not be eliminated if only reduced indices are considered. A cached matrix can be 
reused for all weights, as long as the set of restriction vectors S does not need to 
be enlarged and the precision B does not need to be increased. For an example on 
how to cache these matrices, the reader is referred to the files at [Raul 2 a] which are 
related to Section [HI 

Finally, we consider restrictions to Jacobi forms of higher rank index. It is not 
necessary to restrict the Jacobi forms that we wish to compute to sz, z G C for 
vectors s G Z^. One can equally well restrict them to sz, z G C^' for any sufficiently 
large set of matrices s G Mjv',Ar(Zi), A^' < A^. This can be seen directly by applying 
Lemma 14.61 to the obvious generalization of Proposition 14.81 In order to compute 
Fourier expansions of Jacobi forms of index L[s\, s G MAr/^Ar(Z) one has to, eventually, 
restrict to scalar indices. However, it is possible to break up the problem into smaller 
pieces, which can be treated in parallel. Also, since restriction to Jacobi forms of 
higher index pertains more information, it is conceivable that the working precision 
B can be decrease. For lattices of large size this results in a significant reduction of 
the number of short vectors that have to be enumerated. 

The previous observation gains even more importance in the light of the following 
speculation. In some cases there is a geometric interpretation of embeddings C "-^ C. 
Such an interpretation possibly allows for computing the restriction map on Fourier 
expansions more efficiently than by bare enumeration of vectors in £. This kind of 
tricks could be helpful when computing Jacobi forms for, say, the Leech lattice. 



7.2. Input data to |GKR11) . In |GKR11] . we used the presented implementation 
to compute Fourier expansions of weakly holomorphic vector valued modular forms 
that serve as input data to the multiplicative Borcherds lift. Input data to Borcherds 
products for the orthogonal group C = U ®U ® Co{—l), where U denotes the hyper- 
bolic plane, corresponds to Jacobi forms of positive index £o- This is a relatively easy 
application of the presented algorithm, since Lq is positive definite, and so we do not 
need to apply Algorithm 12.31 Because weakly holomorphic vector valued modular 
forms are connected to vector valued modular forms via multiplication of sufficiently 
high powers of the modular form A, one can compute them without major difficulty. 
At |Raul2a] . we provide a Sage script that computes weakly holomorphic modular 
forms of non-positive weight k and type Pco{~i)- Such forms are uniquely defined by 
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their principle part 

a{m, fi) q"" tf,. 

/iGdisc Co 
0>mGQ 

We illustrate how to use the script. 

load weakvv_process . sage 
k = -1 

L = matrix (2, [2, 1, 1, 2]) 
p_part = { (-1 , (0,0)): 1 } 
B = 5 

vector_valued_niodular_f orni_by_principle_part (k , L, p_part , B) 

The principle part is pass as a dictionary, whose keys have the following (m, /i). The 
first component is a rational number and the second component is the elementary 
representation of an element in disc C. The values of that dictionary are the Fourier 
coefficients a{m,fi), m < 0. The function will fail if any m is non- negative. It will 
also fail, raising ValueError, if the desired vector valued modular form does not exist. 

8. Vector valued Jacob: forms for lattices of rank 1 and 2 

Vector valued elliptic modular forms associated to cyclic discriminant groups are 
of particular interest, since they are the most basic case. Among the cyclic ones, 
the discriminant groups Dm := (Z/2mZ, V H-> /2m), m > are the simplest. 
They come from positive definite lattices (2m). Consequently, the corresponding 
vector valued modular forms of type p£, can be obtained as components of the theta 
decomposition of holomorphic Jacobi forms of index m. An algorithm to compute 
their Fourier expansions has been known for a long time. A first version was given 
in |EZ85] . and a much more refined version, which is also used in Purple Sage [S^lla] , 
was presented in |Sko84] . 

From a structural point of view, the discriminant forms Dm '■= (Z/2mZ, f i— )■ 
f ^ / 2m) , m < are equally simple, but it is harder to compute Fourier coefficients 
of associated vector valued modular forms. One can employ a theta decomposition 
to relate such vector valued modular forms to skew-holomorphic Jacobi forms as 
defined in |Sko90| . Apart from Eisenstein series and theta series, no construction 
of skew-holomorophic forms is known. The major complication when dealing with 
skew-holomorphic Jacobi forms is that the product of two of them is no longer a 
skew-holomorphic Jacobi form. Thus the Fourier expansions of skew-holomorphic 
Jacoobi forms are much harder to compute than those of holomorphic Jacobi forms. 
Algorithm 12.31 and 14.31 provide us with a relatively easy way to compute the Fourier 
coefficients of vector valued modular forms of type pOm- The Jacobi indices that arise 
during these computations have size 7x7. As explained in Section HJ enumeration 
of short vectors in such lattices is time consuming and Sage does not come with an 
optimized implementation. For this reason we have precomputed several cases that 
are possibly interesting for further applications. 

Using hardware provided by the Max Planck Institute in Bonn, we have computed 
Fourier expansions up to precision 20 of vector valued modular forms of type p£,^ for 
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all > m > —4. A short excerpt is shown in Table [T] and El The complete data is 
available at |Raul2a] . We illustrate how to access and use it. Recall that the elements 
of Mk[p£)) have a Fourier expansion of the form 

0<neQ 

The following code only requires Sage |S"'"llb] to be installed. We consider the case 
m = 2, A; = |, and we want to know the value of a/(0, 4), where / is the basis element 
of M5(p£)_2) that we have computed. We simply type 

fs = load ( vector_valued_f ourier_expansion__2k_5__ni_ -2__ . sobj ' ) 
(fs [0]) [0] [4] 

This gives us —312. The code does the following: We load a list f s of basis elements 
for the space M5(pi3_2)- our case this list has length 1. The elements of this 
list are dictionaries whose keys correspond to elements in Z /2|m|Zi, which label the 
components of modular forms of type They are integers < r < |m|. Note 

that the Fourier expansion of any element of M^^pDm) supported on + c_^, 
yU G -D in the notation of Section |3l Hence the available values suffice to reconstruct 
all Fourier coefficients a/(/i,n) for all /i. 

We have also precomputed Fourier expansion of vector valued modular forms for 
negative definite lattices C of rank 2: 

(l 2) ' (02)' (1 4) ' (04)- 

These expansions are related to |GKRll) . where they serve as input data to an 
algorithm computing Borcherds products. The corresponding weakly holomorphic 
vector valued modular forms can be accessed as was described in Section 17.21 For the 
given indices all elements of M_i{pc(-i)) with pol order at infinity less than or equal 
to 5 are provided. 

9. Special divisors on modular varieties of othogonal type 

The data provided in the previous section can be used to deduce linear equiva- 
lences of special divisors on modular varieties of othogonal type. Let £ be a lattice 
of signature {2,n), n > 3. The orthogonal group 0(£) of C consists of all linear 
transformation of £ (g) Q which preserve C. The discriminant kernel 0(£) C 0(£) is 
the kernel of the action of 0(£) on disc C = / C. A modular variety of orthogonal 
type is of the form 0(£)[disc£]\02,n / O2 x 0„. We will give a model for this variety 
that allows us to easily define special divisors. The reader is referred to |Bru02] for 
details on the construction that we explain. Set 

V, = {z e C0C : qc{z) = 0, {z,z)c > O}, 

where 1 means the complex conjugate in the second tensor component. It is well- 
known that T>e has two connected components. Choose one and denote it by "D^. 
The projectivisation V = /(C \ {0}) is a homogenous space for 02,n- We write 
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m 














1 


2 


3 


4 


5 




1 


1 

4 


5 
4 


9 

4 


13 
4 


17 
4 


21 

4 


as ( 


m, u) 


1 


-70 


-120 


-240 


-550 


-528 






-10 


-48 


-250 


-240 


-480 


-480 


09 ( 

2 


m, yU) 


1 


242 


2640 


11040 


30962 


65760 






2 


480 


4322 


13920 


39360 


73920 


" 13 


(m,/i) 


28 





-433680 


-4736160 


-21626280 


-74216064 


2 




-195 


-58344 


-933595 


-6583560 


-30703920 


-96379920 


(2) 


(m,/i) 





56 


240 


-1440 


704 


-960 


2 




-1 


-120 


-9 


1320 


240 


-5040 



Table 1. Fourier coefficients afc(m, /i) of basis elements of the spaces 
of vector valued elliptic modular forms of weight k and type P(-2)- 



Xn for 0(£)[disc The linear extension of (-,■)£ to £ (g) C allows us to define 

divisors in D^, and D as follows. Given A G set 

Z[\) = {zeV : (^,A),: = 0}. 

This is an analytic divisor on T>. For /i G disc£ and < m G Q, the special divisor 
of discriminant (m, n) is defined as 

(9.1) ^(m,/i)= 5^ Z{\). 

The special Picard group PiCsp(X£) of is 

m, |U)] / {div(/) : / meromorphic on X£,div(/) C Z[Z(m, /i)]}. 

Assume that C = U ®hl ® £'(—1), where lA is the hyperbolic plane. In this case 
its is known that |Bru02l IBFlOj . relat ions linear equivalences of Z{m,fi), that is the 
relations in PiCsp(X£), are given by coefficients of vector valued modular forms. More 
precisely, let be a basis of M2+n (pr) with Fourier expansion 

/iGdisc C mGQ 

Given a finite sum Z = Xl/,£disc£ J2o<m&Q fJ')Z{rn, /x), we have ~ Z G PiCsp(X£) 
if and only if 

(9.2) ai^m, fi) b{m, fi) = 

/iGdisc C' 
0>meQ 

for all 1 < 2 < dim(M2+:i(p£)). 
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16 


25 
16 


41 

16 


57 
16 


73 
16 


89 
16 


105 
16 


121 
16 
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1 

4 


5 
4 


9 
4 


13 
4 


17 
4 


21 

4 


25 
4 


29 
4 




3 


1 

16 


17 
16 


33 
16 


49 
16 


65 
16 


81 
16 


97 
16 


113 
16 
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2 


3 


4 


5 


6 


7 






2 


-48 


-144 


-192 


-716 


-384 


-864 


-768 






-25 


-121 


-192 


-336 


-528 


-624 


-864 


-1321 






-10 


-48 


-250 


-240 


-480 


-480 


-1210 


-720 






-1 


-48 


-144 


-337 


-384 


-673 


-816 


-864 









-92 


-96 


-288 


-384 


-672 


-576 


-1152 
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-92 


-96 


-288 


-384 


-672 


-576 


-1152 


2 




25 


121 


192 


336 


528 


624 


864 


1321 






-10 


-48 


-250 


-240 


-480 


-480 


-1210 


-720 






1 


48 


144 


337 


384 


673 


816 


864 






2 


-48 


-144 


-192 


-716 


-384 


-864 


-768 



Table 2. Fourier coefficients ak{m, fi) of basis elements of the spaces 
of vector valued elliptic modular forms of weight k and type P(-8)- 



In the case C = (m), we say that Xc is a paramodular variety of genus 2 and 
level m. The computations in Section [8] can be used to get precise relations of special 
divisors on these Xc- For example, we have 

Corollary 9.1. On the paramodular variety ^w®w®{-8) of index 4 the following rela- 
tions generate all relations that hold for Heegner divisors Z{m,fi) with m < 2. 

-Z{1, 4) = 100 Z{0, 0) + 3 1) + 3 2) + 3 3), 
Z(i|,3) = 2Z(^,l)-2Z(^,3), 

2) = 78 Z(0, 0) + 3 1) + 3 2) + 3 3), 
Z(f|,l) = 5Z(^,l)-4Z(3^,3), 

-Z(l, 0) = 80 Z(0, 0) + 4 1) + Zil 2) + 2 3), 
-Z(0, 4) = 29 Z{0, 0) + 1) + 3 2) + 3 3), 
-5 3) = 30 Z{0, 0) + 1) + 3 Z{\, 2), 
-5Z(i,2) = 50Z(0,0) + 2Z(^,l). 
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Proof of Corollary \9. il We have to consider vector valued modular forms of weight 
|, whose coefficients are given in Section [SI It is easy to check that relation (19. 2p is 
satisfied in all case, and that they span the space of all relations over Z. □ 

More such relations can be computed using the data presented in Section [SI The 
question of how to compute modular forms with divisors as in the above corollaries 
was considered in [GKRll] . 
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